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1 Introduction 



The cosmological equations (without cosmological constant) for a uniform universe 
with the space-time metric 



-dt 2 + R 2 (t) 



dr 2 



kr 2 



+ r 2 {d6 2 + sin 2 



are 



'dR\ 8ttG - , 



r .... d 

R -r = t 
dt dt 



R 3 (pc 2 +p) 



1.21 



;i.s) 



where R = R(t) scales the comoving coordinates, pit) is the energy density, p(t) 
the pressure, k = 0, ±1 the curvature constant and t is a universal cosmic time. 
The Friedmann models |2| have the properties: p ~ 0, k — 0, ±1 (A; = is in fact 
the Einstein-de Sitter model 0), and were improved and made more realistic by 
Lemaitre [Q who included non-interacting radiation. If one denotes by subscripts 
m and r matter and radiation, respectively, one finds from (1.3), for non-interacting 



matter and radiation: p m oc R , p r oc R , for p m << Pr, 
becomes 



a r R + a m R 



Pr 




p r c 2 /3. Eq. (1.2) 



;i.4) 



a r = 87rGp r i? 4 /3, 



a,, 



8nGp m R 3 /3. 



;i.s) 



Eq. (1.4) is just the Lemaitre equation. Its solution for k = 1 was given by de Sitter 
0, Tolman and Alpher and Herman 0, for fc = 0, ±1 by Chernin || and Cohen 
@, and for k = by Jacobs pi. Cosmological models containing both matter and 



radiation were also discussed by Mcintosh |TI] and Harrison \V1 . 

For non-interacting matter and radiation, p r oc T 4 and hence T oc -R -1 ; and p m oc 
n oc T 3 , where n is the mean matter number density. Therefore, the dimensionless 
quantity 

V = n(^X = n^(-4--) 3 (1.6) 



k R T. 



k R T(°) 



is constant. Let us evaluate this quantity by assuming that we are dealing with 
baryonic matter. 

There are about 10 57 nucleons in a typical star. There are about 10 11 galaxies in 
the universe, and each galaxy has about 10 11 stars. Thus there are about 10 79 baryons 
in the universe (in comparison with 10 89 photons, a number which is obtained by 
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thermodynamic arguments). The present size of the observable universe is 10 28 cm. 
The baryon number density is therefore given by 

10 79 

Ub ~ ^(1Q28)3 ~ 10 6 Cm 3 , 

which we take for n<®. Taking ~ 2.7 K (the microwave background), one obtains 
7] ~ 10~ 9 . Let (3 = p^l Pm ( zero superscript denotes the present epoch); then because 
Pr c 2 ~ k B T(k B T/hc) 3 , 

h t(°) 

/?~^-~l(r 3 , (1.7) 
7]m n c z 

where m n is the nucleon mass. Let the superscript (1) denotes the epoch where 
Pr-Pm] then 

— = = (3F6 \ (1.8) 

and p« = p$P~ 3 ~ 10 9 p, ( n } - A1 so, = T(°)//9 ~ 3000 K, and when T > T«, the 
energy density of radiation exceeds that of matter. 
The numerical values of the energy densities are [I3| 

p^ ~ 2 ■ 10~ 31 g cm" 3 , pf ] ~ 4 ■ 10~ 34 g cm" 3 , (1.9) 

pW ~ 10~ 22 g cm" 3 . (1.10) 

It follows that the temperature dependence of the energy densities can consistently 
have the form 

p m (T) ~ 10- 31 Q gem" 3 , p r (T) ~ 10- 22 (3^3) gem- 3 . (1.11) 

As remarked by Harrison [12], when T > ~ m e (? jks — 5 • 10 9 K (m e is the 
electron mass), the Lemaitre equation breaks down because of lepton and hadron pair 
production. At temperatures T ~ the interaction of the hadrons is strong and 
has mainly a resonant character, the masses of the resonances being comparable with 
the temperature. Under these conditions hadronic matter is neither an ideal nor an 
ultra-relativistic gas, and can be well characterized by a resonance spectrum |TJJ] 

r(m) ~ m a exp (m/T ), (1.12) 

where a and T are parameters. The statistical sum diverges when T > T , which 
indicates that the theory involves a limiting temperature (the so-called "hadronic 
boiling point" of Hagedorn [14J]), whose numerical value is found to be 

T = m w c 2 /k B ~ 10 12 K (1.13) 
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(m^ is the pion mass), and in the vicinity of To particle creation is so violent that it 
is impossible to exceed this temperature. In the temperature range ~ T ~ T 
the hadronic gas is thought of as being well described by the "realistic" equation of 
state suggested by Shuryak [JL5|], p, p oc T 6 (which gives the value of the velocity of 
sound c 2 = dp/ dp = 0.20), in agreement with some of the relevant experimental data 
16 , |l7l , |18|] . There are indications that c 2 may, however, take the values both more 



and less than 0.20, i.e., 0.25 [19, 20, 21 



18, 22, 23, 24 



or 0.17 

The existence of the limiting temperature T suggests that there is some phase 
transition. It also means that the Lemaitre equation, although reasonably realistic, 
is a good approximation throughout the lifetime of the universe except for its earliest 
moments. 

It is clear that, in fact, the known resonances form an essentially discrete set 
of states, and the well-known arguments applied above are based on an approxi- 
mate idealization that considers the envelope of these resonances as a continuous 
mass spectrum. There is, however, a consistent (proper time) formulation of a mani- 
festly covariant statistical mechanics |25|, EO, 1271], based on the ideas of Fock p8i and 



Stueckelberg p9f , in which the four components of energy-momentum are considered 
as independent degrees of freedom, permitting fluctuations from the mass shell. 

In the present paper we shall use this manifestly-covariant framework (which 
we review briefly in the next section) as a model for the description of phenomena 
which take place in hot hadronic matter. We argue that the phase transition at 
the Hagedorn limiting temperature can represent a phase transition from an off- 
shell sector of the theory, in which the relativistic ensemble is described within the 
framework of the manifestly covariant relativistic statistical mechanics mentioned 
above, to an ultrarelativistic independent particle phase. Recently we have studied 
thermodynamic properties of the off-shell phase |27, 30[] and its possible consequences 
in hadronic physics PH, astrophysics |52fl, and cosmology |33], p4| . 
As shown in ref. 



the behavior of hadronic matter below the Hagedorn limiting 
temperature coincides with that of a system which includes both particles and an- 
tiparticles, with the additional mass potential |2J| px — 0; such a system, within the 
covariant framework, is described by the equation of state corresponding to Shuryak's 
"realistic" one. In the present paper we shall show that off-shell matter, if included 
into the Friedmann-Lemaitre model on an equal footing with baryonic matter and 
non-interacting radiation, has energy density comparable to that of non-interacting 
radiation at temperatures of the order of the Hagedorn temperature, indicating the 
possibility of a phase transition from strongly interacting (off-shell) phase to non- 
interacting one. 
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2 Relativistic TV-body system 



In the framework of a manifestly covariant relativistic statistical mechanics, the dy- 
namical evolution of a system of N particles, for the classical case, is governed by 
equations of motion that are of the form of Hamilton equations for the motion of N 
events which generate the space-time trajectories (particle world lines) as functions 



of a continuous Poincare-invariant parameter r, called the historical time |29|, |35 
These events are characterized by their positions g M = (t, q) and energy-momenta 
p^ = (E, p) in an 8iV-dimensional phase-space. For the quantum case, the system 
is characterized by the wave function ip T (q 1 ,q 2 , . . . , g^v) € L 2 (R m ), with the measure 
d i qid i q 2 • • ■ d^q N = d AN q (qi = gf ; fj, — 0, 1, 2, 3; i — 1, 2, . . . , N), describing the dis- 
tribution of events, which evolves with a generalized Schrodinger equation The 
collection of events (called "concatenation" j36|) along each world line corresponds 
to a particle, and hence, the evolution of the state of the iV-event system describes, 
a posteriori, the history in space and time of an iV-particle system. 

For a system of N interacting events (and hence, particles) one takes ]35| (we use 
the metric = (-, +, +, +)) 

^ = Ew + y(9l,92 '---' 9jv) ' (2 ' 1} 

i 

where M is a given fixed parameter (an intrinsic property of the particles), with the 
dimension of mass, which we take to be the same for all the particles of the system. 
The Hamilton equations are 

dq- = dK_ = 

dr dp^ M' 

dtf = _0K_ = _ dV_ _ 

dr dq iit 

In the quantum theory, the generalized Schrodinger equation 

d 

l g-Vv(gi, 92, • • • , Qn) = Kip r {qi, q 2 , . . . , qx) (2.3) 
describes the evolution of the iV-body wave function if) T (qi, q 2 , ■ ■ ■ , <27v)- 



2.1 Covariant thermodynamics and cosmology 



Thermodynamic functions for a many-body system can be derived from the grand 
partition function, which, for an ensemble of off-shell events at temperature T is 
defined by the following expression, modified, in comparison with the standard one, 
by the presence of the term (j,rK (henceforth we use the system of units in which 
k — c — Ub — 1, unless the other specified): 



Tr 



cxp 



(2.4) 
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where p is the chemical potential and px is the additional mass potential of the 
ensemble ||25|| . Here E^ is the operator of the total energy of the iV-body system, 
£{N) _ Y,f = iPi, K^ N ' is the generalized A^-body Hamiltonian defined by Eq. (2.1), 
and N is the operator of the number of events (and therefore, particles). As shown 
in ref. |33| (we shall review this point briefly below), the quantity — px K , where k 
is the generalized Hamiltonian density, k = (K^ N ')/V, represents a (55)-component 
of a generalized energy-momentum tensor, which in the local rest frame takes on 
the form T a p = diag(— p,p,p, p, — upx^), a,/3 = 0,1,2,3,5. Here the sign of a 
corresponds to the invariance group of the extended (x^, r) manifold which could be 
0(4, 1) (a = 1) or 0(3,2) (a = —1.) In the present paper we shall restrict ourselves 
to the case a = 1, i.e., the 0(4, 1) invariance of the r) manifold having the metric 
g^= (-+,+,+,+). 

Expressions for p and p, using the grand canonical ensemble obtained by Horwitz, 
Schieve and Piron P5J in their study of manifestly covariant statistical mechanics, 



were found in BO] in terms of confluent hypergeometric functions as 



An A p\ ^ s 2 2p K T 
where Tav is a characteristic interval of r for a trajectory to pass through a small 



representative four-volume. For T small, one finds from (2. 5), (2. 6) 0, |33[ that 
p, p oc T 6 , p ~ 5p, and, in fact, that px^ oc T 7 is negligible in comparison with p. 
On the other hand, for T large, it follows from these expressions that p, p, px^ oc T 5 , 
p ~ l/4p ~ —pKK- Now, we remark that the energy- momentum tensor 

T*" = {p + p)u»u y - pg^ (2.7) 

can be extended to a five-dimensional form 

^T aP = ( ^% u , ^T 55 ) ; (2.8) 

the requirement that the limiting case of the corresponding gravitational theory (for 
zero curvature in the r direction) coincide with the Einstein equations results in the 
identification = ^T^ v and ^T 55 = —p^K. For high temperature, it therefore 

follows that (as discussed in 



T a/3 = {p + p)u a u p - pg a/3 , u x u x = -1, (2.9) 



so that, in the local rest frame, T a p = diag (—p,p,p,p,p)- 
The first law of thermodynamics reads 



dE = TdS - pdV + p l dN i , (2.10) 
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where S is the entropy, V = V(t) (the universal cosmic time t corresponds to the r of 
(2.2) and (2.3)) is a comoving volume element of an ideally uniform universe, and iV, is 
the number of the ith kind of particles having the chemical potential /ij, which implies 
that the energy E is a thermodynamic function of S,V, and iVj : E = E(S,V, iVj). 
In the case of varying iVj it is convenient to consider the set of variables (T, V, //j) 
instead of (S, V, Ni). The change of variables is obtained with the help of a Legendre 
transformation 

n(T,V,(i i )=E-TS-fj^N i , (2.11) 

which introduces the thermodynamic potential Q. The use of (2.10) in the differential 
of (2.11) leads to 

dtt = -SdT - pdV - NidfM. (2.12) 
The coefficients S, p and iVj are given by the partial derivatives 

fdtt\ fdQ\ fdQ\ 

8= -{*>L ; > = -wL ; w «=-y,v- (2 - 13) 

A fundamental result of statistical mechanics relates the thermodynamic potential to 
the grand partition function: 

n(T,V,(n)=-ThiZ. (2.14) 
By combining (2.13) with (2.14), the entropy and the pressure become 

S ^ P.15) 

In the case of no interactions the latter relation reduces to 

p =l\nZ = -^, (2.17) 



which is the equation of state of a free relativistic ensemble |25| 

^ = lnZ. (2.18) 

Substitution of Q — —pV into Eq. (2.11) yields the thermodynamic relation 

E = -pV + TS + fHNi, (2.19) 
which represents the formula for the entropy |37] 



s = E + pV-fMjNi = (p + p-fMjN^V (2 2Q) 
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where p = E/V, Noi = Ni/V are the corresponding energy and particle number 
densities. 

Since the lepton and baryon numbers are relatively small and therefore negligible 
in the early universe, equilibrium conditions correspond to 

fjii = 0. (2.21) 

Then the first law reads 

d{pV) = TdS - pdV, (2.22) 
and Eq. (2.20) takes on the form 

S =^K. (2.23) 
Since V ~ R 3 , it follows from (1.3) and (2.22) that 

and the entropy is conserved. Note that Eq. (1.3) represents the local energy conser- 
vation. Indeed, for the energy-momentum tensor of an ideal cosmological fluid, 

T" v = {p + p)u»u u - pg^, u p u p = -1, 

the local energy conservation, V M T MZ/ = 0, takes on the form 

R tt +Hp + P) ll = > (2 ' 25) 

which reduces to (1.3). 

By rewriting (2.22) in the form 

TdS = d[(p + p)V] -Vdp (2.26) 

and using the integrability condition P7| 

d 2 S d 2 S 
dTdV ~ dVdT' 
the energy density and pressure are related as follows, 



(2.27) 



P + * = T%. (2.28) 

Given an equation of state, the temperature dependence of p and p can be derived with 
the help of (2.28). Assuming that in the early universe the pressure is proportional 



to the energy density and using Zeldovich's equation 

P=(T-l)p, (2.29) 
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we find 

p,pcxT r /( r - 1 ). (2.30) 
In the standard framework one can obtain from (2. 25), (2. 29), 

p,potR~ 3r . (2.31) 

For the radiation-dominated universe T = 4/3, so that 

p,p(xT A (x R~ A . (2.32) 

In the covariant framework Eq. (2.28), and therefore (2.30), hold, while Eq. (2.25) 
should be modified. For the generalized energy-momentum tensor (2.8) the local 
energy conservation, V a T alS = 0, yields 

i?f + 4(p + p)^ = 0, (2.33) 

resulting, via (2.29), in 

p,potR~ Ar . (2.34) 
Since in this case [^3J T = 5/4, we obtain from (2. 30), (2. 34) 

p,p ex. T 5 oc R~ 5 . (2.35) 



3 "Realistic" Friedmann model 

Let us modify the Friedmann-Lemaitre model by inclusion of off-shell matter with 
the equation of state[]p m ' oc T 5 oc -R -5 , p m > = p m 'C 2 /4. The modified equation reads 



'dR\ 2 
dt 1 



a m >R~ 3 + a r R~ 2 + c^iT 1 - k, (3.1) 
a m > =8nG Pm/ R 5 /3, (3.2) 



ni 



a r ,a m being defined in (1.5). 

Let = p,/p^\ where, as previously, the superscript (1) denotes the epoch of 
T« ~ 3000 K^ecause p^c 2 ~ k B T(k B T /Tic) 4 , 

, k B T<n k B TW/{3 * fl T<°> 

Pr V PmC 2 /(3 6 PmC 2 



1 The equation of state of relativistic off-shell matter p = 1/4/3 was obtained also by Hakim 
within a different framework. 
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If superscript (3) denotes the epoch where p m > ~ p r (implying that the pressure of 
radiation and off-shell matter are equal; from (1.11) one sees that the contribution of 
baryonic matter is negligible at high temperature), then 



a. 



RW = P'RV, (3.4) 



and 

$ = pW)" 4 ~ 10 36 P« ~ 10 14 g cm" 3 . (3.5) 

Also 

T(3) = ~ 3 . 10 12 K. (3.6) 

P 

We see that this temperature is just of the order of the Hagedorn one. At T ~T&\ off- 
shell matter and radiation have comparable energy densities, indicating the possibility 
of a phase transition. In general, 

p v (T)^10 14 (^-^) 5 gcm- 3 . (3.7) 

Since the cosmological phase transition at T c ~ 10 12 K is normally associated with 
the transition from a strongly interacting hadronic phase to a weakly interacting 
quark-gluon plasma phase |24|, |40| , we associate the T 4 -phase above the transition 
temperature with a phase of of weakly interacting quarks and gluons, and T 5 -phase 
below the transition temperature with a phase of strongly interacting hadrons (as we 
show in [41 1, T 4 behavior indeed results from the very high temperature asymptotic 



behavior of the distribution function, for small In this phase, particles undergo 
continual mutual interaction and are necessarily off-shell. Therefore, the effect of 
strong interaction in such a system may be represented by the off-shellness of its 
particles, justifying the use of the off-shell framework for the description of this phase. 

If the phase transition at temperature ~ is not sufficiently sharp, there exists 
some region of the "mixed" phase, in which both (or the three, including baryonic 
matter) phases coexist (a similar situation may occur for both hadronic and the 
quark-gluon plasma phases, according to the "realistic" scenario suggested by Shuryak 
||1 5|| ) . If the off-shell phase extended well below T^ 3 \ it would have energy density 
comparable with that of baryonic matter at T ~ 10 8 K, as follows from (1.11) and 
(3.7), so that the temperature range of the mixed phase would be 10 s K ~ T ~ 10 12 K, 
where p m ~ p m > ~ p r . In fact, the mixed phase, as we explain below, may exist in 
the temperature range 

10 10 K ~ T ~ 10 12 K, (3.8) 

where 10 10 K is the electron-positron threshold for hadron production in the e + e~ 
annihilation process, and at temperatures well below 10 10 K the off-shell phase is 
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exponentially suppressed. For T from the range (3.8), we obtain, with the help of the 
relations 

2 dp dp/dT 
dp dp/dT' 

P=Pm+Pr+Pm'' P = Pm + Pr + P m > , 

the following expression for the velocity of sound: 

n2 _ n m p m + n r p r + n m >p m > 
n m p m + n r p r + n m >p m > 
n m c 2 m + n r c 2 r p r /p m + n m > c\ p m > / p m 

-, (d.yj 



c 



n m + n r p r /p m + n m 'P m '/Pv 



where n m , n r , n m ' are the powers of temperature in the corresponding formulas for 
the pressure and the energy densities, and (? m = dp m /dp m = p m /Pm, etc. are the 
sound velocities in the corresponding phases. Using the relations 



we finally obtain 



C m 1) C r g' C m ' ^: 



2 4/3 p r /pm + 5/4 p m >/p r , 



C 2 ~ ' ^ ' "l^m rrn _ ( 3 .1 ) 

3 + Ap r /p m + 5p m '/p m 

If the energy densities of the corresponding phases were of the same order, p m ~ p r ~ 
p m /, one would obtain 

31 

c 2 ~ -0.21. (3.11) 
144 v ' 

In fact, the use of Eqs. (1.11), (3. 7) in the formula (3.10) provides for the sound 
velocity a numerical value which is practically 1/3, up to temperature ~ 2 • 10 12 
K. At T = (~ 3 ■ 10 12 K), Eq. (3.10) yields c 2 w 0.28. One sees that for the 
model of the mixed phase, a dip in the sound velocity as a function of temperature 
in the vicinity of the transition temperature is predicted, in agreement with both the 



phenomenological models [|T5|, [42f and the recent data on QCD lattice simulations 



f43| . Let us also note that the value of the sound velocity c 2 = 0.28 at T ~ 10 12 K 



was reported in ref. pH 



The value of the sound velocity c ~ 0.25, in agreement with the relevant exper- 



imental data |L9], |2y, ^TJ, is obtained only in the assumption of a sufficiently sharp 
transition, when the phase below the transition temperature is that of strongly inter- 
acting (off-shell) matter alone, with temperature dependence T 5 . 

Thus, the "realistic" Friedmann model studied here describes qualitatively the 
decrease of the sound velocity in relativistic gas at values of T in the vicinity of 10 12 
K, which reaches the numerical value c 2 ~ 0.25, consistent with at least some of 



the experimental results [19, 20, 21 1, and significantly less than the ultra-relativistic 
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Stefan-Boltzmann gas value 0.33. It may be thought that the progressive population 
of states of higher mass creates an admixture of heavy particles, thus decreasing the 
factor relating pressure to energy density. 

The above estimates have been made in a classical framework. Now we shall 
show how similar results can be obtained within the framework of the conventional 
(on-shell) quantum theory. As mentioned above, at ~ m e ~ 5 ■ 10 9 K, lepton 
and hadron pair production starts; in the temperature range 

T(2) < t ~ T°, where 

T is the temperature of the Hagedorn phase transition, T ~ T^ 3 \ the system is 
characterized by a resonance spectrum of the Hagedorn form (1.12), for which there 
is some experimental basis |T4] , |15|j . At T ~ 10 10 — 10 11 K << T , one can neglect the 
exponential in (1.12) and use, as done by Shuryak |T^, fj5| , a resonance spectrum of 
the form 

r(m) ~ m a , (3.12) 

where one chooses on phenomenological grounds a = 1 [ |15|| . 

The expressions for the pressure and energy density are then written as |15, f|5| (we 



neglect the difference in properties of Bose and Fermi particles at high temperature) 

p = 9 - [°°dmr{m) f - = - e -^^/ T = g [°°dm r(m)p m , (3.13) 

3 J a y ' J (2tt) 3 VA; 2 + m 2 J a V ; V ; 

d 3 k 



p = g dmT{m) / -— ^-Vh? + m 2 e~ Vk2+m ' 2/T = g / dmr{m)p m , (3.14) 
Ja J (2n) ,i Ja 

where g is the number of hadronic degrees of freedom, p m and p m are expressions 
for the pressure and energy density of relativistic gas of particles with given mass m, 
and we may choose, sufficient for our present purposes, A = 2m e ~ 10 10 K. Since 
at temperatures T ~ A one should take into account electron-positron annihilation, 



leading, via quark loops, to the formation of hadronic jets ffi> 



e + e — > qq — » hadrons, (3.15) 

which can be thought of as resonances on the excitation background of virtual leptonic 
pair states, the e + e~ threshold should be taken as the lowest possible mass for the 
resonance spectrum r(m), justifying the introduction of the corresponding cut-off in 
the formulas (3.13), (3. 14). 



The use of the standard expression p7 



m 2 T 2 T ^ fm\ 
Pm = ^rK 2 (-) (3.16) 



2tt 2 VT 

in Eq. (3.13), in which we take r(m) = Cm, C = const, yields 
gCT 2 3t ^ fm\ gCT 3 



f . (3. 17 ) 




12 



At temperatures T ~ 10 11 K where T >> A (but still T << To, so that the use of 

justify 

K u (z) ^ - I (/')(- I . :«l 



(3.12) instead of (1.12) is justified), one uses the asymptotic formula fl49 | 

1 



2 v ' V2 



and obtains 



P 



AgC 



71' 



Similarly, one can obtain for such temperatures E3| 

20#C^ 6 



7T^ 



bp, 



so that 



0.20. 



(3.18) 
(3.19) 

(3.20) 
(3.21) 



Equations (3.19)-(3.21) represents the "realistic" equation of state suggested by Shuryak 
for hot hadronic matter ||15|| . 



For T << A, one uses another asymptotic formula |ff9 

K v {z) ~ J ^~eT z , z»l 

V Zi % 

and obtains 



(3.22) 



V 



(3.23) 



gCT^A^ _ A/T 

2 5/2 n 3/2 e 

and similar expression for p. Thus, at temperatures well below A ~ 10 10 K, the 
off-shell phase is suppressed by the exponential. 

The use of r(m) ~ m a in the formulas (3. 13), (3. 14) will analogously lead to the 
equation of state at T ~ 10 11 K [ 15|, 45 1 

p,p~T a+5 , p = p/(a + 4), (3.24) 

which coincides with (3.19)-(3.21) for a — 1. The values for the sound velocity in the 
hadronic phase 0.25 and 0.17, in agreement with some of the experimental results 
18|]-|24|, is achieved, in view of (3.24), for a = and a = 2, respectively. 
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4 Solutions to cosmological equation 

The solutions of Eq. (3.1) for t > t (3) , R > R {3) are well known fll2l 
k = : 



R 

t 

k=l: 
R 

t 

= -1 : 
R 



y/a^T + -a m r , 

1 >— 2 1 3 

-y/a^r + —a m T ; 



/ • • 2 T 

'ot r sm r + a m sin — , 



2 r 1 

2-y/ovsin — | — a m (r — sinr) 



2 2 



^air sinh r + a m sinh 
r 1 



t = 2-y/a^sinh — + -a m (sinh r — r). 



2 2 

For < t < t( 3 \ < R < RW, since the curvature term is negligible in the early 
universe, we have the equation of the general form 

2 

| =a n R 2 ~ n , a n = 87iG P R n /3, 



dR 
~dt 



with n — 5. Eq. (4.1) has the solution 



R = ai^ n 



V(n-2) ( 1! ± T \ 



so that 



and 



—a. 
n 



R 



l/(n-2) / ra ~^ r ) 



2/n 



a 



l/n+2/r 



t = - 

n 



2 /SttG A" 1 / 2 



(4.1) 

(4.2) 
(4.3) 

(4.4) 
(4.5) 



For the energy density at T = p, ~ 10 14 g cm 3 , we obtain from (4.5), with 
n — 5, 

t {3) ~ 10" 5 c. (4.6) 



Equation (4.1) corresponds to that of Zeldovich and Novikov [50| for the early charged 
symmetric universe, whereas the Lemaitre equation (1.4) is for the subsequent charge 
asymmetric universe in which the asymmetry is either local or global. 
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5 Concluding remarks 



We have studied the cosmological Friedmann model modified by introduction of off- 
shell matter with the equation of state p oc T 5 oc -R -5 , p — p/4. The energy density of 
such a matter is comparable with that of the non-interacting radiation at temperature 
of the order of the Hagedorn limiting temperature, ~ 10 12 K, indicating the possibility 
of a phase transition. 

A cosmological phase transition at T c ~ 10 12 K is normally associated with the 
transition from a strongly interacting hadronic phase to a weakly interacting quark- 
gluon plasma phase |24|, f40fl . The simplest classical model considered in the present 
paper implies the possibility of a first order phase transition from the T 5 -phase of 
strongly interacting (off-shell) matter to the T 4 -phase of non-interacting radiation- 
like matter. Cosmological consequences of such a phase transition are discussed in 
ref . |34j] , where it is suggested that the transition may be sufficiently smooth (second 
order) to preserve the expansion rate. Although a first order phase transition might 
be preferable for some cosmological implications, due to the fluctuations which are 
generated at the transition and could produce planetary mass black holes |51| which, 
in turn, could provide a possible explanation for the dark matter of the universe and 
even be seeds in galaxy formation |52], [53| , experimental indications on the order of this 
phase transition are still absent. Indeed, presently available lattice data on SU(N) 
pure gauge theory lattice simulations indicate that a phase transition to a weakly 
interacting phase is of apparently first order for SU (3) and second order for SU (2) 
theory [[54]]. In ref. |55| , however, the apparent first order nature of the transition 



in the case of SU (3) pure gauge theory has been called in question. Moreover, there 
are indications from lattice QCD calculations that when fermions are included, the 
phase transition may be of second or higher order | 56fl . In this remarked by 

Ornik and Weiner |[24|| , the phase transition would be hardly distinguishable from a 
situation in which no phase transition would have taken place (radiation-dominated 
universe alone). 
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